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Abstract 

The low-lying energy excitations of 2-flavour QCD in the chiral limit confined 
to a small spatial box ((5-regime) are that of an 0(4) rotator. In this work, we 
calculate the corrections due to the presence of a nonzero quark mass up to NNL 
order by means of dimensional regularised chiral perturbation theory. The final 
result for the energy gap of the system only involves the low-energy constants 
F, Ai, A2 and B. 



1 Introduction and summary 



Quantum chromodynamics, the quantum field theory describing the interactions 
between quarks and gluons, is non-perturbative at low energies (large distances). 
An expansion in the coupling constant breaks down, since the latter becomes 
large (of 0{1)) at low energies. In order to study QCD properties at low energies, 
we have to apply other methods. 

The lattice provides the only known non-perturbative treatment of QCD in 
terms of the fundamental fields, the quarks and the gluons. However, lattice 
QCD calculations are very time-consuming, and so we are restricted to relatively 
small volumes or relatively large lattice spacings. Thus, the results we obtain 
from lattice QCD simulations suffer from finite size effects and discretisation 
errors. 

In this work, we will restrict ourselves to two light quark flavours in the 
isospin limit (m„ = = m). QCD is known to be spontaneously broken 
in case the u and d quarks are massless. The QCD Lagrangian is invariant 
under the global SU{2)l x SU{2)ii transformations, the vacuum of QCD is 
only invariant under global SU{2)v transformations. According to Goldstone's 
theorem [T] , the breaking of a global, continuous symmetry gives rise to massless 
particles, the Goldstone bosons (GB). 

At low energies, QCD can be described by chiral perturbation theory (ChPT) 
[H El S] 7 an effective field theory with pions as degrees of freedom. The low- 
energy coefficients (LECs), associated to the operators, showing up in effective 
Lagrangian, have to be determined by experiments or by measurements from 
lattice simulations. Since ChPT can also be applied if the system is enclosed in 
a finite volume ^ |6j , finite volume effects can be calculated analytically. 

We consider a volume V — Ls Lt, where the temporal extent 

of the box and the Compton wavelength of the pion are both much larger than 
the spatial size of the box, Lt ^ Lg, MLg ^ 1 respectively. This special regime 
is called 5-regime (^-expansion) [7J. Here, M refers to the leading term for the 
pion mass in infinite volume 

= 2mB . (1) 

Since SU{2)l x SU{2)]i is the covering group of 0(4), we use a 0(4) nonlinear 
sigma model in d — 3-1-1 Euclidean space-time dimensions to describe the 
effective theory. 

Due to the large Compton wavelength of the pion, large compared to the 
spatial extent of the box, the system exhibits a global mode which slowly rotates 
in the internal space. These slow modes and the related energy excitations 
interest us in this work. In the chiral limit (m — ^ 0), the low- lying energy 
excitations of this system can be described by the spectrum of the 0(4) rotator 

Ei=^-^, ^ = 0,1,..., (2) 

where this form is only valid up to NNL order. Here, £ denotes the 0(4) "angular 
momentum" (isospin) and O the moment of inertia which is given as = F^L^ 
at leading order [7, 8, 9J. The low energy constant F is the pion decay constant in 
the chiral limit. Allowing for higher order terms in the effective Lagrangian the 
moment of inertia will receive fast mode corrections proportional to \/[FLs)^ 
at NL order [lOj, respectively proportional to l/(FLs)^ at NNL order [TT]. We 
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will denote the dimensionless expansion parameter as 

1 



(3) 



throughout this work. 

For sufficiently small quark masses M < l/F'^L^, the system is still domi- 
nated by the rotator spectrum. The latter then obtains small corrections due to 
the symmetry breaking terms. The dimensionless expansion parameter related 
to the symmetry corrections will be denoted by 

r = F^LIM . (4) 

The NLO correction to the rotator spectrum coming from the symmetry break- 
ing terms is proportional to r'* [7] . 

In Section [31 we will illustrate how to separate the slow modes from the fast 
modes. In order to apply perturbation theory, the effective action has to be 
expanded in terms of the fast modes. A detailed discussion for the symmetry- 
breaking terms is given in Sectional In Section [5l we will integrate out the fast 
modes. The calculations are performed in the zero temperature limit {Lt — ?> cxd), 
and we use dimensional regularisation (DR). This setup is in accordance with 
calculations in 

The remaining quantum mechanical 0(4) rotator, now in a external "mag- 
netic" field, is characterised by the moment of inertia, and by the symmetry- 
breaking, external field ?], where Q and rj contain corrections from the fast 
modes. 

Finally, we apply perturbation theory in quantum mechanics in order to 
determine the corrections to the unperturbed rotator spectrum. Since the two 
expansion parameters enter at NLO as (5^ ^ and r"* respectively, we consider these 
two combinations to be of the same order ~ = 0{6^). We are interested 
in corrections up to NNLO, therefore we have to take into account any terms 
up to 0{d'^). 

We quote the final result for the energy gap E^^ up to NNL order 
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15 



120 152 



(5) 



where is given in [11] as 



1 



2G* 



0.088431628 



5o'G*^(ilog(AiL,)2 +log(A2L.)2 



and 



77 = f^lIm"^ ( 1 - 



■iG* 
2FHA 



G* and OqG* are given by the following numbers 

G* = -0.2257849591 
d^G* = -0.8375369106 . 



(6) 



(7) 

(8) 
(9) 



The low-energy constants li and I2 enter ^ over their intrinsic scales Ai, re- 
spectively A2 112]. 
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2 The effective action up to 0{p'^) 



The low-energy properties of QCD with two nearly massless flavours, considered 
in infinite volume, can be covered by an effective field theory in terms of the 
pseudo-Goldstone bosons, the pions. We choose to describe the GB dynamics 
by an 0(4) nonlinear sigma model. The effective Lagrangian is ordered in a sys- 
tematic way by the increasing number of derivatives and increasing power of the 
symmetry-breaking field H. Derivatives are counted as 0{p) and the symmetry- 
breaking field is counted as 0{p^). Using the convention which identifies the 
term A(^d,h) in the effective action containing operators with d derivatives and 
the field H to the power h, we write 

Acs — -4(2,0) + "4(04) + -4(4,0) + "4(2,1) + -4(o,2) + • ■ • , 
for the effective action up to 0{p'^). The explicit expressions for these terms are 

^2,0)= Jdx ^d^S{x)d,Six) , (10) 

A^a^i-)^- jdxY,HS{x), (11) 

^(4,0) = Jdx ) (d^S{x)d^S{x)) ' + igf ) (^dj{x)dj{x)) \ (12) 

A24)= Jdxh^(^HS{x)) (d^S{x)d^S{x)) , (13) 

A^o,2)^^ Jdxk2^[HSix)y . (14) 

S{x) is a 4-component vector of unit length, S{x)'^ = 1, and we choose the 
external field H to point in the zeroth direction H — {H, 0, 0). Terms which 
do not depend on x at all have been omitted, since they will enter as an overall 
factor in the path integral representation. We have used the same conventions 
as in [IB] for the low-energy constants. 

Due to the explicit symmetry breaking, the pions acquire a mass. If we 
consider the external field H to be small, the pseudo-Goldstone boson mass is 
given by 

= (15) 

at leading order in H. In the following, we will replace the combination "EH 
by M^F^. By comparing Eq. with Eq. ([1]), we can identify the external 
symmetry breaking parameter H with the quark mass m and E with BF^. 

The effective action p^ - ([CT) is based on conventions from condensed matter 
physics. In chiral perturbation theory, the pion fields are usually parametrised 
by 5/7(2) matrices, and the low-energy constants in the 0{p'^) effective La- 
grangian are labelled by I2, Is-, and hi. We determine the relations between 
the LECs by comparing the 0{p'^) effective action in fT31 with the terms PU)) - 

m 

= -4/1 , ) = -4/2 , 

ki =Ia, k2 — h + h , ks = hi . 
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3 Separating the fast modes from the slow modes 



Chiral perturbation theory can also be applied in finite volume [5J [HI. In the 
(5-regime {Lt ^ Ls and MLs <C 1), the fields S'(a;)on a given time slice are 
strongly correlated and exhibit a net "magnetisation" 



rn(t) = — jdx S{t, x) , 



rnit) = m{t)e{t) . 



(16) 



The direction e{t) of the net "magnetisation" performs a slow rotation in the 
internal space. These slow modes have to be treated non-perturbatively. The 
fluctuations (fast modes) around direction of the "magnetisation" can be in- 
tegrated out in perturbation theory. Therefore, the fast modes have to be 
separated from the slow modes. 

We incorporate the collective behaviour of the variables S by introducing 



1 = Yl Jdmit)S^^^ 



m{t) 



1 



dx S{t, x) 



into the partition function of the system. The partition function then reads as 



Z =Y[ ldS{x)6 
Y[ fdrn{t)6^ 



S^x) - 1 



m(t) 



1 



dx S{t, x) 



(17) 



By choosing appropriate field redefinitions, which have been worked out in |11| . 
Z can be written as 



Y ldxTv{t,x) 



It is indicated that in the effective action the variable S{x) will be replaced 
by the combination 

S{t, x) = n{t)V^{t)R{t, x) . (19) 



The 4-component vector R of unit length is parametrised as 
R{x) = (v/l-7r2(a;),7r(x)) . 



(20) 

Eqs. ([20| and (fT8| ensure that the slow modes are not a part of the 7r-fields. In 
fact the k — (/cq, fc = 0) modes have to be left out when we calculate the Green's 
functions. Furthermore, it follows that 

y Jdxn,{t,x)=0, i = l,...,N-l. (21) 

ri(t) and V{t) are 0{N) matricefl The first column of the matrix n(t) is 

e^{t) = n^o{t) , a = 0,...,A^~l, (22) 



^In paper the notation was introduced to our V{t). We decided for this change, 
since the quark condensate is denoted, here (and in many other works), by S. 
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and the matrix V{t) has the fohowing structure 



V{t) = 





: V{t) 

V / 



(23) 



The partition function (|18p is expressed as a path integral over the slow 
modes e{t) and as a path integral over the fast modes tt{x). The small fluc- 
tuations (fast modes) can be integrated out in perturbation theory. Therefore, 
we have to expand the effective action in terms of the 7r-fields, up to the de- 
sired order. The 7r-fields enter by replacing the variables S according to the 
parametrisations defined in Eq. and Eq. (PO)) . 

We have not defined the two matrices n{t) and V{t) completely. For the 
symmetry- breaking terms the unknown (undetermined) parts of fl(t) and V(t) 
will drop out. Hence, the fixing of these two matrices will not be discussed in 
this work. A detailed treatment of this issue is given in 

4 The symmetry-breaking terms up to NNLO 



The terms pl))) - (ri4l) of the effective action are still expressed in the field variable 
S. As mentioned before, we have to expand the effective action in Eq. (|18p in 
the TT-ficlds up to the desired order (NNL). In this paper, the focus is set only 
on the symmetry-breaking terms, i.e. the terms ^(o.i), -^(2,1) a-nd .4(0,2) • The 
expansion of the symmetric terms has been covered in (IT] . 

4.1 ^(0,1) up to NNLO 

The symmetry-breaking term of the 0{p^) effective action is given by 



A(^o.i)iS)^~F'M' JdxSoix). (24) 

Using the relations ^ to express S in terms of 7r-fields, So 

can be written as 

So{x) ^ noc,{t)Vpc{t)Rp{x) 

= eo{t)RQ{x) + no^{t)Vj,{t)^Tj{x) 

We use the convention that Greek indices run from to — 1, whereas Latin 
indices run from 1 to A^ — 1, and we will stick to this convention throughout 
this paper, as long as nothing else is mentioned. 

The term proportional to '!Tj{x) contains some unknown elements of the 
matrices V and fl. According to (|2ip this term will, however, vanish when Sq 
is plugged into Eq. ([2^ . After having expanded Ra{x) in terms of the 7r-fields, 
.4(0,1) reads as 



A(^OA) = -F'M^ Jdxeoit) (^l-^7v^{x) + ..}j . 



(25) 
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The dots indicate that there are terms of higher order in the 7r-fields. However, 
we can truncate the expansion already at this order, since the leading contribu- 
tion from the symmetry breaking will enter already as a NLO correction in the 
rotator spectrum. 



4.2 ^(0,2) and ^(2,1) up to NNLO 

The symmetry-breaking terms entering at 0{p'^) are given by Eq. (|14p . and 
Eq. respectively. We will give arguments why contributions from these two 
terms will only enter beyond NNLO. 

In order to get a crude estimate for the sizes of the corrections coming from 
these two terms, we consider the system in a simplified form. We neglect the 
fast mode corrections. The variable S{x) can then be replaced by e{t), and the 
leading contribution from the symmetric term simplifies to 

^2t3 ^ ' 1 



^(2,0) ~ F'Li jdt-e{t)e{t) (26) 

which can be interpreted as the action of the 0(4) rotator. The energy eigen- 
values of this system are proportional to 1/ F^L^. 

Under the same assumptions (tt — > 0), the symmetry breaking term .4,(0,2) 
is reduced to the following form 

« " Jdtk^M^L^eoitf , (27) 

This term can now be regarded as a perturbation to the symmetric rotator (M 
is small). The field eo is quantity of 0{1), naively the spectrum of the rotator 
will receive a correction proportional to ^ M'^L^ 

^ + 0{M^Ll) = (1 + OiM^LlF')) . 

The relative correction to the unperturbed rotator energy is then proportional to 
F'^L^M^. The latter combination, however, is of 0{d^). This can be shown eas- 
ily by expressing the strength of the perturbation in terms of the dimensionless 
expansion parameters S and r 

F^LlM^ = r^je = oiS"") . 

Similar considerations lead to the same conclusion for the term .4(2.1)- Ne- 
glecting any fast modes corrections, the corresponding part of the action reads 
as 

dt ^^k{t)k{t)2h— eo{t) . (28) 

This expression can be considered again as a small perturbation to the un- 
perturbed rotator spectrum. Assuming the field eo, which is of 0(1), to be a 
constant, the action reduces to the symmetric case. The corrections to the mo- 
ment of inertia is then proportional to The latter combination is of 0(6'''), 
and therefore also these corrections are beyond NNLO. In addition, we have 
completely neglected the fact that this term would not contribute in the first 
order of the expansion, due to odd number of fields involved. Thus, the actual 
corrections would be even smaller than 0{S'^). 
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5 Integrating out the fast modes 



We consider again the partition function (|18p and insert the effective action, 
which has been expanded in terms of the fast modes, in the exponent exp(— ^cff ). 
While keeping the kinetic terms for e{t) and tt{x) in the exponent, the remaining 
terms are expanded in a Taylor series and will be treated as interactions 



N — l 



• exp 



mm 



■ exp 



F^ 



2 ^tj V 



(t) 7r,(x)7rj(a;) + .. 



(29a) 
(29b) 



F 



^TT'^{x)doTri{x) - (n{x)doTT{x)'^Tri{x) 



1 



^(Tv{x)d^,Tv{x)^ (n{x)df,n{x)^ + 



(29c) 
(29d) 



1 - 



-,(2) 



1 - 94 



(3) 



^e{t)e{t)df,n{x)df,n{x)+ 
+ \Qio{t)Q joi.t)doni{x)donj{x) 



l-e{t)e{t) doTv{x)don{x) + \Qio{t)Qjo{t) 



(29e) 



1 

• (j)oTTi{x)do'Kj{x) + ^ dkTriix)dk'n-j{x)] 
1 + I F^M^eoit) 



fc=i 



(29f) 



l--7r2(x) + 



+ 



2! 



F^M^eo{t)eo{t') 



l-l7v\x)-^7v'iy) 



(29g 



X y 



At this point the fast modes can be integrated out in a systematic way, keeping 
only terms up to 0{S'^). In Eq. (|^^ we used the convention x = {t,x), and 
y ~ {t',y) respectively. Furthermore, it is assumed that after having integrated 
out the fast modes, the limit e has to be taken. The symmetric terms (|29ap - 
(|29fl) have been discussed in [llj, where also the origin of the e is explained in 
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more detail. In the following, the focus is set on the additional symmetry- 
breaking part ( |29gp . 

The pairing of two 7r-field components which live at the same space time 
point X will give a contribution of the form 

(^.(a:V,(x)}cx||Z?*(0). 

D*(0) denotes the finite volume propagator evaluated at x = 0. The * indicates 
that the k = (fco, fc = 0) modes have to be left out when calculating the prop- 
agator. In = 4 dimensions, the finite volume propagator is proportional to 
L^^ . The pairing two 7r-fields results in a contribution 5^ 

= -^G* , (30) 

where G*, given in Eq. ([5|), is the finite volume Green's function evaluated at 
a; = in DR. 

Let us consider expression ( |29g| separately, that means we ignore cross terms 
with any of the symmetric interactions for now. Again, we integrate out the 
fast modes, and we keep only terms up to 0{S'^). Hence, the partition function 
can be written as 



+ 1 jdt |dt'(^^^L^Af2)^(^l-3^)eo(t)eo(0 + ...|. 



(31) 



where, up to NNL order Q is given as quoted in Eq. © The expression inside 
the curly brackets in Eq. (PT|) . can be written as an exponential again 



Z oc Jl y de(t) exp - j'^t ^^(^)^(*) - Veo{t) 



(32) 



Here, we introduced the new expression 77 which we have already given in Eq. ([7]) 
as 

3G* 



?7 = F^LlM^i 1 - 



Since the leading contribution from the symmetry-breaking interactions is a 
0(15^) correction to the rotator spectrum, the corrections to 77 are only consid- 
ered up to (FLs)~^. 

A similar argument can be used to neglect any cross terms between the 
symmetric interactions and the symmetry-breaking interactions (on the level 
of the fast modes). The only possible cross terms which would connect fast 
modes from the symmetric part with fast modes from the symmetry-breaking 
part will enter only beyond 0{5'^). Thus, the cross terms factorise in a simple 
way and we can write the symmetric part and the symmetry-breaking part as 
an exponential ([5^ . 

The original problem, described by an effective field theory in d = 4 Eu- 
clidean space-time dimensions, has been reduced to a 1-dimensional system 
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described by the action 



(33) 



where the variables e(t) satisfy the constraint 



e{t)e{t) = 1 . 



In the chiral hmit (77 — 0), we identify the system ([55)1 as that of a sym- 
metric 0(4) quantum mechanical rotator, where up to NNLO, the spectrum 
is given by In the explicitly broken case, the rotator is considered in an 
external "magnetic" field which points along the zeroth direction, r/ is assumed 
to be small, so that the symmetry-breaking potential can be treated as a per- 
turbation. The corrections to the unperturbed rotator are then calculated in 
perturbation theory in quantum mechanics. Considering again corrections up 
to 0{S^), requires perturbation theory in quantum mechanics up to 4th order. 

6 The quantum mechanical 0(4) rotator 

In this section, we will do the discussion for general > 3, although we are 
explicitly interested in case iV = 4, finally. The system described by Eq. ((33)) 
can be identified as an 0{N) rotator in a small external "magnetic" field. The 
Hamilton operator for this constraint system is given by 



where L can be considered as the "angular momentum" operator in the A^- 
dimensional internal space. Again, is the moment of inertia and A, given 



is the small, dimensionless parameter which denotes the strength of the external 
"magnetic" field. 

The energy gap of the system is defined as the difference between the energy 
of the first excited state and the ground state energy. In the chiral limit (A 0), 
the ground state energy is zero and the energy gap is simply given by the energy 
of first excited state. Due to the presence of the "magnetic" field, the energy 
levels of the first excited state split up. For = 4 the first excited state splits 
up into a singlet and a triplet. The triplet provides the lower energy difference 
and is identified as the energy gap. Up to 0{X^) the energy gap has already 
been calculated in [7], but without taking into account fast modes corrections. 

Besides taking into account also fast modes corrections, we are interested in 
corrections up to 0{X^). Since A is assumed to be small, we write the energy 
spectrum of as a power series in A up to 0{X'^). 




(34) 



by 



A = r,e , 



(35) 




(36) 
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Here, e^^l denotes the spectrum in the chiral hmit (A 0) which for general 
A'^ > 3 is given by 

42 =4°^ = ^^(^ + ^-2), ^ = 0,1,.... (37) 

In order to determine the coefficients s^^], we apply standard perturbation theory 
in quantum mechanics. 

The energy eigenstates of the unperturbed system can be represented by 
associated Jacobi polynomials 

|£,fc)^p(f (z). (38) 
The polynomials -P/^'' (z) are defined as 

'N-l\ / ri \ ( 



(39) 



{z), (40) 

and have been normalised to Pf ^\l) = 1. k denotes the 0{N — 1) "angular 
momentum" and runs from = 0, . . . ,i. The variable z takes values in the interval 
[—1,1] and can be identified with the direction of the "magnetic" field in the 
internal A''-dimensional space. For fixed N and k, two of these polynomials 
arc orthogonal in the interval [—1,1] with respect to the weighting function 
«;(z) = (1 - z2)(^-3)/2 

1 

= Jdz (1 - z^)'^PS\^)P'e^^(^) - • (41) 
-1 

In perturbation theory matrix elements of the following form have to be 
calculated repeatedly 

, . (£,fc|z|f,fc) (zP/.f (z)Pf,\z)) 

where only matrix elements of the form Vf^f^-^ for fc = 0, . . . , ^, and Vf^g_^ for 
k = 0, . . . ,i — 1 are nonzero. As a consequence, the corrections to the orders 
0(A) and 0{\^) vanish 

4!i=45=0, (43) 

and at a given order of A^", n — 1,2,..., only a finite number of terms survive 
the (in general) infinite sums for the coefficients e^^-*, e^^-* respectively 

T/fe -irk T/fe T/fe 

(2) ^ ^e-it , ^i+u / . .X 

^.fe ^(0) _ ^(0) ^ JO) _ (0) ' 



\rk -t/k T/fe T/fe T/fe T/fe T/fe T/fe 

^''^ ~ (JO) _ (0) X 2^.(0) _ ^(0) N + / (0) _ (0) .2 (JO) _ JO) s 



r/k T/fe T/fe T/fe 



(JO) _ JO) X2 (JO) _ (0) X2 



(45) 
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Using Eqs. ([57|) - (|^5)l . the energies E^^k (EH) can be calculated for arbitrary 
TV > 3. 



7 The energy gap for = 4 



The energy gap for = 4 is defined as the difference between the energy of first 
excited state = 1, fc = 0, 1) (singlet, triplet) and the ground state {£ — 0,k — 
0). In Tab. [TJ we quote the coefficients for the ground state energy correction 
and for the first excited state energy correction. The triplet {£ — l,k — 1) 
provides the lower energy difference, and therefore the energy gap is defined as 



El, = Ei^i{\) — £'o,o(A) . 
Up to O(A^), the explicit result for energy gap reads atH 



1 A^ 193 
^ 15 " 120 152 



(46) 



(47) 



Inserting the results for 8 and 77, according to Eq. and Eq. ^ respec- 
tively, we recover the final result in Eq. ([S]). The leading term in A is given by 
the dimensionless expansion parameter r^, and up to 0{5'^) A reads as 



A^ ^ F^L^Hl^ 



1 



7G* 



(48) 





N 


4 








^i,k 







k = 








1 

6 


5 
432 


e = 


1 


k = 





3 
2 


1 

15 


317 
27000 






k ^ 


1 


3 
2 


1 

15 


23 
27000 



Table 1: Here, the coefficients e^";, (n = 0,2,4) are given for the ground state {£ = 
k — 0), for the singlet {(. — l,k — 0), and for the triplet {£ — k = 1). 



8 The constraints on and M 

The formula for the energy gap Eq. ([5]) involves two different expansions. The 
requirement that theses two expansion are applicable will put some constraints 
on the values of Lg and M. These values can be considered as a rough estimate 
on the domain of Lg and M, where the approximation used above is valid. 

^The NL contribution is in agreement with [7], as long as we neglect the fast mode correc- 
tions, i.e. El, = 3/{2F'^Ll) [l + {F^Ll^M'^)/15 + ...]. 
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^NNL 


Mi[MeV] 


M2[MeV] 


M2Ls 


2.0 


0.59 


-0.10 


99 


213 


2.16 


2.5 


0.38 


-0.05 


80 


109 


1.38 


3.0 


0.26 


-0.03 


66 


63 


0.96 



Table 2: This table shows the NLO {9nl) and the NNLO (Onnl) corrections to the 
moment of inertia for some selected values of Ls . Mo is simply Mi = 1/Ls, and M2 is 
defined according to H51|) . As numerical input we used F — 86.2 MeV, Ai — 120 MeV 
and A2 = 1200 MeV, obtained from jTS]. 



What is the constraint on Lg? The moment of inertia O receives corrections 
proportional to (5^ at NLO, and proportional ot S"^ at NNLO respectively 

e = F^Ll [1 + 0NL + 0nnl] . (49) 

In order to have this expansion working properly, the corrections 0nl and ^nnl 
should be small. In fact, we require that 6'nl should be roughly 50% or smaller. 
This requirement puts a lower bound on L^. An estimate for the lower bound 
of the box size can be obtained from Tab. [51 There, the size of the corrections 
at NLO and at NNLO are given for some selected values of Lg- For Ls — 2.0 fni 
the NLO corrections turn out to be rather large, about 60%. Going to towards 
Ls = 2.5 fm the NLO corrections decrease to 40%. Thus, we conclude that 
volumes about Ls > 2.5 fm should be considered. 



What are the constraints on Ml The (5-regime requires the mass to be 
much smaller than the inverse box size. We define Mi to be 

Ml = ^ . (50) 

This relation gives a first, upper bound on the mass for a given Ls- The corre- 
sponding values are quoted in Tab. [21 Taking Ls = 2.5 fm as a reference value 
for the lower bound on Ls, the mass should be smaller than M — 80 MeV. 

However, this constraint does not account for the fact that the expansion for 
the energy gap ([S]) can brake down. This issue can be circumvented requiring 
the first correction to the energy gap (r'*/15) to be small, e.g. 50% or smaller. 
We define M2 as the second upper bound satisfying the relation 

15 2' ^ ^ 

for a given Ls- It turns out that this constraint only becomes relevant for larger 
values of Ls, i.e. around Ls = 3fm. 

From the considerations above we conclude that an appropriate choice for 
Ls and M is crucial, in order to be in the domain where the two expansions 
work properly, and the formula for the energy gap is valid. We defined the lower 
bound on Ls to be about 2.5 fm. This forces the leading order term in the pion 
mass M to be smaller than 80 MeV, and therefore, we have to use quarks with 
masses below the physical quark masses. Since M scales roughly as ^/L^, the 
upper bound for the mass decreases even faster with increasing Ls- 
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